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Abstract 

We study the weakly interacting Hubbard model on the square lattice using 

a one-loop renormalization group approach. The transition temperature T c 

between the metallic and (nearly) ordered states is found. In the parquet 

regime, (T c 3> the dominant correlations at temperatures below T c are 

antiferromagnetic while in the BCS regime (T c <C \fj,\) at T c the d-wave singlet 

pairing susceptibility is most divergent. 
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A theoretical understanding of systems with two different kinds of coexisting and possibly 
competing instabilities of density-wave (DW) and superconducting (SC) type remains one 
of the central problems of the theory of interacting fermions. In particular, the question 
of a purely electronic mechanism of superconductivity, induced by an incipient instability 
of the density-wave type, remains difficult if one wants to go beyond the most qualitative 
level. Typically the existence and strength of the DW correlations and their coupling to 
superconducting pairing depend strongly on the dimensionality and on the geometry of the 
Fermi Surface. 

Experimentally, superconductivity in the vicinity of an insulating and/or spin density 
wave state is a well-known property of the phase diagrams for several families of fermion 
systems: for example, in the quasi-one-dimensional Bechgaard salts, Jl| superconductivity 
(possibly of d-type 0) replaces an SDW state as one increases pressure, due to two ef- 
fects: decreasing of the nesting properties of the Fermi surface, and suppression of umklapp 
scattering between electrons. The phase diagram of the quasi-two-dimensional organic su- 
perconductors of the ET family also shows a pressure induced SDW-SC transition ||, but 
the symmetry of the SC order parameter is not yet clear. || In the high-T c superconduc- 
tors a few percent doping transforms an insulating antiferromagnet into a superconductor 
f|, probably of the d x 2_ y 2 symmetry. || In two dimensions the simplest model showing an 
insulating antiferromagnetic state already at weak coupling is the Hubbard model, often con- 
sidered as a "minimal model" for high-T c superconductors. Although extensively studied, 
the question of doping induced superconductivity in the vicinity of the antiferromagnetic 
state in this model still remains unanswered. In the weak coupling limit the problem of the 
interdependence of different kinds of instabilities can be treated using renormalization group 
methods [§,f|, as has been successfully done for the case of quasi-one-dimensional systems 

raj- 

In the present paper we present a one-loop renormalization group analysis of the Hub- 
bard model. Because of its perturbative nature this gives a quantitatively correct phase 
diagram only for weak coupling and does not allow us to make precise predictions about 
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realistic systems where the coupling is rather strong. However, we are able to describe the 
physical mechanisms which lead to different instabilities of the model. Provided that there 
is no transition to a qualitatively different regime for strong repulsion, this should give at 
least some qualitative insight even outside the strictly perturbative regime, at least for the 
effective low energy model ||. Recently we have formulated one-loop renormalizat ion-group 
equations for the Hubbard model, |§ based on the requirement that the vertex (the effective 
interaction) must be invariant under changes of the energy cutoff A about the Fermi sur- 
face. This way of renormalizat ion is known as field theory approach. || We distinguished two 
regimes, separated by a crossover energy T co which is a function of the chemical potential fi. 
In the parquet regime, A > T co , the contributions to the renormalizat ion from the particle- 
particle and particle-hole loops are both important. In the BCS regime, A < T co , only the 
particle-particle loop contributes to the flow, while the particle-hole part is negligibly small. 
In the parquet regime both loops behave like log 2 A, while in the BCS regime, the p-p loop 
crosses over to log A and the p-h loop disappears with some positive power of A. For the 
exactly half-filled case, the parquet equations gave an instability in the antiferromagnetic 
channel []12[ . The same result was found within a simple scaling theory [ITS] which takes 



into account only processes between electrons at the van Hove points in the corners of the 
square Fermi surface. Introducing the chemical potential as a cutoff for the p-h term of the 
flow equation this theory also gives a transition to d-wave superconductivity. 

An important aspect of the renormalization group is that in many cases A can be inter- 
preted as the temperature to logarithmic precision. This approximation consists in renor- 
malizing the vertex at some temperature T only by virtual processes involving "quantum" 
electrons, those with energy larger than T. Then the vertex depends only on energy variables 
greater than T and its dependence on energies inside the area ±A about the Fermi energy 
is irrelevant. In that sense the field theoretic approach can give correct thermodynamics 
only for logarithmic-scale-invariant problems, because in that case the internal fermionic 
propagators A, A', B, B', C, and C in fig. |IJa) are all exactly at the energy ±A and, con- 
sequently, the vertex at the step A of renormalization contains only contributions from the 



degrees of freedom outside the ring ±A. If there is no logarithmic scale invariance, as for 
the case when A > T co ^ 0, in the field theory approach the lines A,B and C are on the 
shell, but A', B' and C can even lie inside the ring ±A, and A can thus not be considered 
as an effective temperature. In order to obtain thermodynamics from the renormalization 
group we thus have to modify the bookkeeping of mode elimination in a way so that the 
modes inside the temperature ring ±T about the Fermi energy never enter the integration. 
This can be done using the Kadanoff- Wilson mode elimination technique |14| as given by 



Polchinski's equations. fl5Hl7| m the BCS regime, the field theory approach is valid because 
there only the p-p channel survives and scales just as log A. 

In this letter we solve Polchinski's flow equation at the one-loop level, as shown schemat- 
ically in fig.|]. This allows us to find the renormalization of the interaction in the nontrivial 
case when the proximity of half-filling, via nesting and van Hove singularities, makes both 
DW and SC tendencies strong and the critical temperature T c can be already in the parquet 
regime. The renormalization of different correlation functions gives us the phase diagram. 

In the Kadanoff-Wilson scheme, the one-loop renormalization flow of the interaction 
U(Ki, K 2 , K 3 , K 4: (K 1 , K 2 , K 3 )) as a function of the three energy-momenta K; t = (cjj,kj) is 
defined graphically by fig.|IJ(a), where the propagators A, B, and C are on the shell ±A, but 
A', B', and C are constrained to run only over states with |£| > A. The interactions in the 
loops are also to be taken at the same cutoff as the propagators A', B' and C and not at 
the actual cutoff A. This means that the renormalization of the interaction is nonlocal in 
A: this is the price we have to pay to get thermodynamics from the renormalization group. 
For the Hubbard model the initial condition is £7a=a = U, where A is the initial cutoff 
which we take to be equal to the bandwidth (= 8t) in order to take degrees of freedom from 
the whole Brillouin zone into account. It is convenient to introduce the logarithmic scale 
I = In8t/A. 

The fact that U\ is a function of six variables makes the renormalization equations very 
difficult. Power counting arguments can provide some drastic simplifications, in a way that 
allows us to eliminate from U\ all irrelevant variables, namely the frequencies uu{ and energy 
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FIG. 1. (a) The one-loop renormalization of the vertex, (b) The Polchinski equation describes 
the exact flow of all vertices. The contractions labeled with I are integrated over a phase space shell 
of energy width Ad A. (c) The 6-leg diagram generated by the first term in Polchinski's equation. 
The second term applied to it generates the loops for the renormalization of the 4-leg vertex, (d) 
The organization of the mode elimination for a non-half-filled case. Dashed lines are the lines of 
constant "angle" 6. 
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variables £k ; of the interacting particles. A problem with scaling in the parquet regime with 
\i 7^ is that if one chooses processes exactly at the Fermi surface as marginally relevant, 
the nesting and umklapp scattering processes appear formally as irrelevant, whereas this 
irrelevancy only really occurs when A i.e. in the BCS regime. On the other hand, 

the irrelevance of the energy variables allows us to choose as marginally relevant a set of 
processes between electrons at any distance from the Fermi surface in the ring ±A, not 
necessarily at the Fermi surface. Thus, in the parquet regime we keep nesting and umklapp 
processes even for fi ^ by considering processes between electrons exactly at the square 
Fermi surface of the half-filled case. Up to irrelevant corrections U then depends only on 
projections of the momenta on this square, as shown in Fig. 1(d), i.e. = Ux{0i, 9 2 , 3 ). 
On the other side of the crossover, in the BCS regime there remains only the BCS amplitude 
Va(0x,9 3 ) = U^{9\ = 9 2 + it,9 3 = 9^ + 7r), but now the particles are at the Fermi surface, 
and not on the square, since the square is outside the ring ±A. 

The renormalization equation for the interaction as function of 0j only is of the same 
form as in ref. |J, but now with 



P ee {U,U}(9 1 ,6 2 ,6 3 ) 



e{u, u}{e x , 2 , e 3 ) + e{xu, xu}(e 1} e 2 , e 3 ) 



(i) 



and 



Mu 1 ,u 2 }{e 1 ,e 2 ,e 3 ) 



n{[/x, U 2 }(6i, 9 2 , 9 3 ) + U{U 1} U 2 }(9 3 , 4 , 00, 



(2) 



where 
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-qj - A) 
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* u i.\^^\(W>h)U2,\^^\{Oi,0,O a ). (4) 

k u is the momentum of a particle at angle 9 with energy £ = uA. J~(e, 9) = J[(x, y)/(e, 9)} is 
the Jacobian of the transformation from rectangular coordinates in momentum space to the 
polar ones. X is the exchange operator defined as XU(9i, 92, #3) = U($2, 9i, 63). U\ and £7 2 
in the expression (^) symbolize U or XU as given in ref. 0. Note that the interactions U, 
Ui, and U 2 in the expressions (|3|) and ((D are not at the actual cutoff A, but at a cutoff given 
by the configuration of external legs and the integration momentum h u (entering via total 
momentum q ee = k x + k 2 and momentum transfer q e ^ = k x — k 3 ). To solve the renormal- 
ization equations we have discretized the ^-dependence of the interaction, approximating in 
that way the function U\(9i, 9 2 , 63) by a set of coupling constants Z7a(*i> h)- shows 
the flow of some of the 93 coupling constants for the case of a discretization in 16 patches. 
There is one pole at l c = ln(8£/A c ), determined by the initial interaction and the chemical 
potential, and in this case it is in the parquet regime. All coupling constants which have 
logarithmic corrections diverge for I = l c . We identify this point as the critical temperature 
T c = 8texp(— l c ) in a mean field sense. The critical temperature at half-filling, T c ° is slightly 
inferior to the mean-field one, but the ratio l^ IF ' fl c = 0.985 does not depend on the strength 
of the interaction. As one sees in fig.§, T c is finite for any filling, in contrast to mean-field 
calculations where T c falls to zero at a threshold doping given approximately by = T c °. 

To determine the dominant fluctuations at T = T c , we calculate the correlation matrices 
for antiferromagnetism x SDW {@ii $2) and for superconductivity X (^i> ^2)1 defined as 

4(01,025 In -r 2 |) = Jje, jje 2 J(e 1 ,9 1 )J(e 2 ,9 2 ) 

(A^(e 1 ^ 1 ;r 1 )A*(e 2 ,e 2 ;r 2 )}, (5) 

with 5 = SC, SDW . The symbols ">" mean that the energy integrals are over energies 
outside of the shell ±A. Consequently, x SC an d X SDW have the interpretation of the sus- 
ceptibilities at the temperature T = A. At the beginning of the renormalization when 
A = A = 8t both susceptibilities are zero. 
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FIG. 2. The flow of a few typical (among 93) scattering amplitudes for a Fermi surface covered 
by 16 patches, for chemical potential [i = 8texp(— 7.8) and initial interaction U = 4t/3. 

The order parameter variables are 

Af (e, 9; t) = J2 ^k^^k+^r), (6) 

Af w (e } 6;t) = Y: ^k(T)*- CTlk+(7ri7r ) +q (T), (7) 

a 

where k is given by the angle 9 and the energy e. We consider only the static and long- 
wavelength limit and follow the renormalization of the maximal eigenvalues of both cor- 
relation matrices; in the SDW channel the corresponding eigenvector belongs to the A\ 
representation (s-wave), while superconductivity has B\ symmetry (d^.^-wave). For low 
doping (small \i) the divergence at A = T c occurs in the magnetic channel (the "SDW" 
region in fig. |3[). For higher doping, the pairing susceptibility diverges first, i.e. the "d- 
wave" region is a d x 2_ y 2 superconductor. The triple point metal-SDW-superconductor is 
at the crossover line T w As we are considering a two-dimensional system here one 
should be careful about the interpretation of T c : in the case of magnetism, this indicates 
the onset of well-defined finite-range correlations. For weak interactions, this is typically a 
very well-defined crossover. |TS|] In the case of pairing T c can be identified with the onset of 
quasi-long-range order. 
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FIG. 3. The phase diagram. 



In our calculations we have neglected self-energy corrections which are in Polchinski's 
formalism given by Hartree-Fock like terms with renormalized u and q dependent vertices. 
The broadening and redistribution of the spectral weight of the quasiparticles is then deter- 
mined by the dynamics of the vertex, which is irrelevant in lowest order by power counting 
and is therefore neglected. One should however notice that at the two-loop order self-energy 
effects do become important, as known from the one-dimensional case. |T(J In that sense, 
our T c should be understood as a temperature where the effects of interaction start to change 
strongly not only the two-particles correlations, but the single particle properties as well. 

Using the above results and arguments, a more detailed description of the phase diagram, 
fig.^, can be given: the transition between Fermi liquid and the magnetic phase occurs in the 
parquet regime. Here the precursor effects, divergence of the commensurate SDW correlation 
function and relevance of umklapp processes, suggest that T c is the temperature of a metal- 
insulator transition occurring together or very close to a magnetic instability. For a more 
specific description of this phase one should include self-energy terms. In the pairing phase 
(d-wave in fig.|3|) one expects a spin-gap, due to the formation of d-wave singlets, as one 
sees from the renormalization the pairing correlations. The maximum pairing temperature 
is about T°/4, in fact quite close to conserving approximation calculations |I9| and to the 
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experimental value of about 1/6 in lanthanum high-T c compounds. || The d x 2_ y 2 symmetry 
of the pairing is also in agreement with the majority of experiments in cuprates. || 

From large N arguments ||,|9| we know that the self-energy corrections will disappear as 
T c /t if T c is deeply in the BCS regime. Consequently, far from half filling ||, T c is a very 
good approximation for a superconducting transition, stabilized already by an infinitesimal 
inter-plane hopping. Finally, mean-field arguments |l7j suggest that one expects an incom- 
mensurate SDW only in the BCS regime, where perfect nesting is impossible. However, the 
precision of our calculation (we cut the Brillouin zone into up to 24 ^-patches) is not sufficient 
to check whether a magnetic correlation function diverges at T c at some incommensurate 
wave vector. 

In conclusion, we have investigated the two-dimensional Hubbard model using a pertur- 
bative renormalization group approach. This allows us to treat instabilities in the particle- 
particle and particle-hole channels on an equal footing and in particular to address the 
important question of superconductivity induced by spin fluctuations. In the vicinity of 
half-filling we do indeed find a sizeable superconducting transition temperature, about one 
fourth of the typical magnetic energy scale. The pairing is of the d x 2_ y 2 type. 
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